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Abstract 

The convergence of no-core shell model (NCSM) calculations using renormalization 
group evolved low-momentum two-nucleon interactions is studied for light nuclei up 
to ''Li. Because no additional transformation was used in applying the NCSM frame- 
work, the energy calculations satisfy the variational principle for a given Hamilto- 
nian. Dramatic improvements in convergence are found as the cutoffs are lowered. 
The renormalization group equations are truncated at two-body interactions, so the 
evolution is only approximately unitary and converged energies for A ^ 3 vary with 
the cutoff. This approximation is systematic, however, and for useful cutoff ranges 
the energy variation is comparable to natural-size truncation errors inherent from 
the initial chiral effective field theory potential. 



1 Introduction 



The nonperturbative nature of inter-nucleon interactions is strongly scale or 
resolution dependent and can be radically altered by using renormalization 
group (RG) methods to decouple low- and high-momentum modes in nu- 
clear Hamiltonians [T|l2l[3lH|5|6|7] . A consequence of this decoupling is that 
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few- and many-body calculations become more tractable at lower resolu- 
tions [8]|9lll0llllp2] , which implies that ab-initio basis expansion methods such 
as the No-Core Shell Model (NCSM) should exhibit much improved conver- 
gence properties. Moreover, the RG approach has the advantage of being able 
to vary the cutoff as a tool to probe the quality of the many-body approxima- 
tion and to provide estimates of the size of omitted terms in the Hamiltonian 
(e.g., three-nucleon interactions) [8]. 

In this paper, we present NCSM calculations of ground states and some low- 
lying natural-parity excited states for light nuclei (^H, ^He, ^He, ^Li and ^Li) 
using RG-evolved nucleon-nucleon (NN) interactions as input. These poten- 
tials include similarity renormalization group (SRG) interactions and both 
sharp [3] and smooth [5j cutoff versions of low-momentum interactions Viowfc 
derived from one of the N^LO effective field theory (EFT) two-nucleon po- 
tentials p^fll] . The dependence of convergence properties on the cutoff is of 
particular interestQin applying the NCSM framework, no additional transfor- 
mation of the Hamiltonian is performed. Thus, the energy calculations here 
satisfy the variational principle for the NN interaction at each cutoff. De- 
tails of the NCSM approach can be found in Refs. P^16lll7|18|19j and pre- 
vious applications of the NCSM to soft bare NN potentials are described in 
Refs. [201I2T1I22] . 

All of the calculations are for NN interactions only. Since the physics is quali- 
tatively incomplete without including at least three-nucleon (3N) interactions, 
these should be viewed as baseline calculations that set the stage for inclusion 
of 3N interactions in the near future [23]; some areas of investigation have been 
deferred until then. The additional computational costs of 3N interactions will 
severely limit the feasible basis sizes for all but the lightest nuclei. So rather 
than use large spaces to ensure convergence of all the NN-only calculations, 
we have used relatively small spaces for all nuclei. With the improved conver- 
gence rates at lower momentum cutoffs, we have been able to perform all of 
the calculations on a small computer cluster. This is sufficient to establish that 
the scaling of computer resources should be favorable for adding three-body 
interactions. 

The RG equations inevitably shift strength from two-body to many-body in- 
teractions as A or A is lowered, even if many-body interactions are initially 
zero. If we only keep the NN part, the transformations are only approximately 
unitary and few-body binding energies will change ("run") with A or A. This 
approximation is not a problem because the initial chiral EFT potentials have 
truncation errors and associated theoretical error bands for observables. We 



The SRG potentials are characterized by a flow parameter A that sets the scale of 
decoupling of low-energy and high-energy physics, thus playing the role of a cutoff. 
For convenience we will refer to both the SRG A and the Viowfc A as cutoffs. 
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can ensure that the variation of energies (and other observables) are compara- 
ble to the EFT truncation errors and that the hierarchy of many-body forces 
is preserved by restricting the range of RG cutoffs. By tracking the running 
of the energies, we can assess the expected net contributions from many-body 
forces in a more complete future calculation. 

We stress at the outset that although the converged (or extrapolated) binding 
energies for certain cutoffs are quite close to experimental binding energies, 
this should not be taken as motivation for using those cutoffs and neglecting 
three-body interactions. Since nuclear matter calculated with low-momentum 
interactions does not saturate with NN only [1], it is clear that the importance 
of 3N contributions will increase as one moves higher in the table of nuclides. 

The outline of the paper is as follows. In Section [2] we look at the conver- 
gence of ground-state energies of light nuclei in a harmonic-oscillator basis 
using SRG-evolved momentum-space interactions. We find that decoupling 
documented for two-body systems in Ref. [24J is also present in heavier nuclei 
(with NN only). The SRG convergence patterns are compared to those for 
V\owk potentials using regulators with varying degrees of smoothness in Sec- 
tion [31 A simple extrapolation procedure, which is needed for the larger cutoffs 
in nuclei with A > 3 and for all cutoffs in ^Li, yields energies with error bars 
that exhibit the net evolution ("running") of three-body (and higher-body) 
contributions. The error bars reflect the NCSM convergence and the cutoff 
variation assesses the impact of missing many-body interactions. These re- 
sults are exhibited in Section HI Convergence for selected excited states and 
radii are shown in Section [51 We conclude in Section [61 with a summary and 
outlook for future calculations. 



2 Convergence of SRG Ground-State Energies 



In this section, we present results for ground-state energies of light nuclei 
for a range of harmonic-oscillator basis sizes and oscillator parameters using 
potentials evolved with the similarity renormalization group [25ll26f27] . The 
SRG flow is a series of unitary transformations parametrized by a flow variable 
s: 

H, = Uis)HU\s)=T + Vs, (1) 

where H is the original Hamiltonian (corresponding to s = 0) and T is the 
kinetic energy, which is taken to be independent of s (thereby defining Vg). 
The choice of transformations leads to the flow equation for Vs, 

^ = [[Tn, Hg], H,] = [[Td, Vs],H,] , (2) 
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which is evaluated here for NN interactions in a partial-wave momentum ba- 
sis Pl24f28j . The operator is a diagonal matrix in this basis. The results of 
Refs. ^f2^ and all results here use Tc = T, but we have also checked conver- 
gence for = and = T^- For = T in momentum representation, 
A = 1/ s^^^ is a more useful flow variable that can be thought of as a cutoff on 
momentum transfers, and we use it exclusively from now on. 

We use primarily the N^LO potential from Ref. [13], which has been the sub- 
ject of previous investigations in the NCSM [TSfTU] . In those investigations, 
Lee-Suzuki transformations were used to derive effective interactions that im- 
proved convergence, but which lose the variational property when making a 
cluster approximation. In contrast, the present calculations use interactions 
evolved in free space, so they simply represent alternative Hamiltonians; thus 
the calculations retain the variational principle. Because the evolution to lower 
momenta yields nearly universal interactions for A < 2fm^^ [DEISIEI, conver- 
gence properties starting from other initial N^LO [H] or phenomenological 
potentials [29] are quantitatively similar. For simplicity, in these first NCSM 
calculations we evolve only the neutron-proton part of nuclear forces, with 
Coulomb added after the RG evolution. 

The choice of To = T in Eq. ([2]) leads to a partial diagonalization of the 
potentials in momentum representation, which is the source of decoupling 
of low- and high-momentum contributions. We expect reduced short-range 
correlations and improved convergence in almost any basis as a result of this 
decoupling. The evolution modifies only the short-distance part of operators; 
for the Hamiltonian this is important to maintain the hierarchy of many-body 
forces in the initial EFT Hamiltonian (as long as we don't evolve too far). 
However, one could also consider tailoring the evolution to a particular basis 
such as harmonic oscillators. In this case, choosing to be the one-body 
harmonic oscillator Hamiltonian (with a two-body term added to the potential 
as usual with the NCSM [T^16|17j ) will yield a flow to more diagonal and 
therefore more decoupled potentials in this basis, accelerating convergence. 
The downside is that there may be serious negative consequences for long- 
distance operators and many-body forces, analogous to what is observed for 
Lee-Suzuki transformations in the NCSM [30J. These issues will be explored 
in a future investigation. 

Matrix elements of the momentum-space potential are evaluated in a harmonic- 
oscillator basis for input to the Many-Fermion Dynamics (MFD) code [31 j . 
which constructs the many-body matrix elements and performs the diagonal- 
ization. The basis is specified by the oscillator parameter Ml and by A^max- 
-^max defines the maximum number of oscillator quanta (increments based on 
sums of single particle quanta, 2n + I) allowed above the lowest many-body 
configuration for a given nucleus. Thus iVmax = 2 implies the inclusion of 
single-particle excitations up to two major oscillator shells or two particles 
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excited simultaneously by one major shell above the lowest-energy oscillator 
configuration. 

We view the present NN-only calculations as setting the baseline for upcoming 
calculations that will include 3N forces. Because the latter will have signifi- 
cant computational costs that will limit the maximum basis size that can be 
used, we study convergence for different A values with at most iVmax = 12 and 
dimension sizes of order 10^ or less. This enables a significant number of ex- 
ploratory cases to be easily and quickly examined on a small computer cluster, 
including the evaluation of a suite of experimental observables for each case. 
The basis dimension for ranges from 24 for A'^max = 2 to 4200 for A^max = 12, 
for ^Li from 800 for iVmax = 2 to 9,692,634 for iVmax = 10 and for ^Li from 
1961 for A^max = 2 to 6,150,449 for A^max = 8. Extending these calculations to 
more powerful computers for including 3N forces is straightforward. 

In Figs. [1] through [6l the ground-state energies of ^H, ^H, ^He, ^He, ^Li, and 
^Li are plotted against the harmonic-oscillator parameter hfl for a range of 
A^max values. The figures with A ^ 4 show results for the initial chiral N^LO 
potential and three SRG-evolved potentials with A = 3.0, 2.0, and 1.5 fm~^. 
For A = 6 and A = 7, our largest spaces were very far from convergence with 
the initial potential, so we show only SRG-evolved potentials with A = 3.0, 
2.0, 1.5, and 1.0 fm'^ 

It may appear in Fig. [T] that the deuteron is converged to different ground- 
state energies for different A values. However, the results simply reflect the 
slow convergence of the weakly bound deuteron in a harmonic oscillator basis; 
all converge with sufficiently large spaces to the same energy as the unevolved 
potential. Only for A = 1.5 fm~^ is the energy converged at the lOkeV level 
by A^max = 12. (For the other cutoffs, there is an irregular convergence trend 
for which A^max = 12 is very close to A^max = 10, but then A^max = 14 will be 
lowered further.) The convergence rate for '^H in Fig. [2] is similar when viewed 
on the same scale. In contrast to the deuteron, the convergence is to noticeably 
different energies in large spaces (indicated by the dotted lines) because the 
SRG evolution truncated to NN-only is only approximately unitary for A ^ 3. 
The spread of energies is a measure of this approximation; for the range of A's 
shown, it is the same order as the truncation error from omitting 3N forces 
in the original Hamiltonian. The running of energies and its implications are 
discussed further in Section HI 

In all cases we flnd rapidly improving convergence with lowered A down to 
A = 1.0 fm~^, which is consistent with an expected decrease in correlations 
in ground-state wave functions. Power-counting arguments [8] imply that the 
hierarchy of many-body forces in heavier nuclei could break down by this point 
(i.e., the 3N and higher-body contributions will become comparable to the NN 
contribution). However, the scaling of many-body forces in the SRG has not 
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Fig. 1. Ground-state energy of the deuteron as a function of hil at four different 
values of A (oo, 3, 2, 1.5 fm^^) for A'^max values up to 12. The initial potential is 
the 500 MeV N^LO NN-only potential from Ref. [H]. The dotted hues show fully 
converged energies obtained with a sufficiently large basis size. 
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Fig. 2. Ground-state energy of the triton as a function of hQ at four different values of 
A (oo, 3, 2, 1.5 fm"^). The initial potential is the 500 MeV N^LO NN-only potential 
from Ref. [13]. The legend from Fig. [1] applies here. The dotted lines show fully 
converged energies from independent A'^max = 48 calculations using a code from 
Ref. M- 
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Fig. 3. Ground-state energy of '^He as a function of /lil at four different values of A 
(oo, 3, 2, 1.5 fm"^). The initial potential is the 500 MeV N^LO NN-only potential 
from Ref. [13]. The legend from Fig. [1] applies here. 
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Fig. 4. Ground-state energy of ^He as a function of HQ, at four different values of 
A (3, 2, 1.5, Ifm"^). The initial potential is the 500 MeV N^LO NN-only potential 
from Ref. 
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Fig. 5. Ground-state energy of ^Li as a function of fi-il at four different values of 
A (3, 2, 1.5, Ifm"^). The initial potential is the 500 MeV N^LO NN-only potential 
from Ref. [H. 
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Fig. 6. Ground-state energy of ''Li as a function of HQ, at three different values of A 
(2, 1.5, Ifm"^). The initial potential is the 500 MeV N^LO NN-only potential from 
Ref. tl3j. 
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yet been established, so the practical lower limit in A for which many-body 
forces are under control is not yet known. Also, it may be different for light 
and heavy nuclei. Once we have established the technology to consistently run 
many-body forces with the SRG (Ref. [23] is a first step in this direction), we 
will have greater freedom to exploit the choice of A in practical calculations. 

In general, the curves of energy vs. HQ and the trends with increasing A^max 
at fixed hfl are more systematic with lower A. For example, compare ^Li at 
A = 3.0 fm~^ to A = 1.5 or 2.0 fm~^. It is also clear that the larger nuclei have 
more regular behavior at any A. This implies that extrapolations to iVmax = oo 
will be more robust for heavier nuclei. 

The trends with smaller A's include increasingly reasonable estimates from 
calculations with smaller A^max values. For example, A^max = 2(4) for A = 
1.5 fm~^ is within 1(0.3) MeV of the converged binding energy for ^He and 
within 5(1.5) MeV for ^Li. To set the scale for how good these estimates are, 
we note that the starting N^LO potential in ^He is at most bound by 1 MeV 
for A'max = 2 and only 8 MeV for A'^ax = 4. For A = 3.0fm~\ A^^^ax = 2(4) is 
still about 14.5(7.3) MeV short of the converged binding energy in ^He. 

At the same time, with smaller A's the larger A^max values converge entirely 
(at the 10 keV level or better) for the lighter nuclei. As A^max increases, the hfl 
dependence gets flatter in a very smooth and systematic way. We also note that 
for a given Amax, the minimum in HQ moves toward lower values as A decreases, 
as is expected for softer interactions. The energies for A = 1.0 fm~^ converge 
remarkably fast. The binding in this case is smaller than for the other A values, 
which is consistent with the running of the net three-body contribution in 
Ref. [8] for a sharp Viovjk interaction in and ^He. Explicit results for the 
running of the net many-body contribution are given in Section |H 

We can compare the convergence patterns observed here with those found in 
previous NCSM investigations. The convergence of Lee-Suzuki effective inter- 
action results with increasing A^^ax is generally non-monotonic at fixed hfl. 
This is seen in numerous examples in the literature (see, e.g., Ref. [22]). It is 
therefore more challenging to make extrapolations in those cases. However, ex- 
trapolations with bare but soft NN interactions, such as JISP16 (Ref. 122]), are 
more straightforward ^33j and similar to the SRG convergence. Extrapolations 
of SRG results are tested in Section HI 

In Ref. [21], the SRG was used to demonstrate how running to lower momen- 
tum decouples low- and high-momentum contributions to matrix elements of 
low-energy observables. Only two-body observables were considered. Using the 
NCSM, we can extend the tests of decoupling with NN interactions to few- 
body systems. One test is to apply a cutoff function to the SRG potential, in 
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the form 

V,{k, k'; = e-^''/'?-^^^"V,{k, k')e~^''''''-^^r , (3) 

where the choice for n controls the smoothness of the cutoff, and then to 
calculate the ground-state energy as fcmax is varied from to cxd. This second 
cutoff (here with n = 8) is simply imposed by hand after the RG cutoff is run 
as a tool to test decoupling: if there is decoupling of low- and high-momentum 
contributions, then we should be able to set matrix elements of the potential 
to zero in a smooth way above some momentum /cmax and still get the same 
answer for low-energy observables. The effect of the second cutoff is shown 
for the triton in Fig. [7] starting from the N^LO potential of Ref. pJi|. For 
comparison, we show in Fig. [H] the same calculation but starting with the 
Argonne fis potential [29] . 

The convergence with increasing fcmax of the binding energy for the bare (A = 
00) interaction to the asymptotic value is determined by the intrinsic cutoffs in 
the potential and reflects the associated decoupling. Thus for AV18, the energy 
is converged to good accuracy by fcmax ~ 7fm~^ while for the N^LO potential 
convergence is reached between 4 and 4.5 fm~^. (Note that the latter result is 
much higher than the naive estimate of 2.5 fm~ based on the 500 MeV cutoff, 
due to a significant high-momentum tail and the associated tensor strength.) 
As the potential is evolved to lower A, the convergence scale for the binding 
energy is set by A. That is, the most rapid changes happen for /cmax up to A and 
then there is a slower approach to the asymptotic value of the energy, which is 
very well converged in all cases by about 1.4 A. These details are common for 
any potential evolved with T to define the SRG [31] . Similar results for other 
nuclei and a quantitative perturbative analysis of decoupling will be given in 
Ref. [31]. 

We have also examined the convergence of the ground-state energy for SRG 
potentials using = and = in Eq. ([2]) rather than the kinetic 
energy T. The convergence behavior in both cases is almost indistinguishable 
from those for T. 



3 Convergence of V\o^k Ground-State Energies 

In this section we compare the convergence for the SRG to that of V\o^k inter- 
actions with sharp and various smooth regulators. In contrast to the SRG 
flow equations for V^, the V\o^k potentials are based on the invariance of 
the low-energy NN T matrix under changes of a cutoff A on relative mo- 
menta [2|3f5] . This is achieved either through coupled differential equations 
for the momentum-space matrix elements of V^ow k or in integral form through a 
free-space Lee-Suzuki transformation [1] . The construction and characteristics 
of Viowfc interactions with smooth regulators is described in Ref. [S]. 
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Fig. 7. Decoupling through the SRG. Binding energy as a function of cutoff parame- 
ter A;max for three different values of A. The cutoff function is exp[— (/c^/fc^g^x)"] with 
n = 8. The initial potential is the 500 MeV N^LO NN-only potential from Ref. [T^ . 
These results include evolved electromagnetic interactions and charge dependences. 




Fig. 8. Decoupling through the SRG. Ground-state energy as a function of cutoff pa- 
rameter /cmax for three different values of A. The cutoff function is exp[— (/c^/fc^g^x)"] 
with n = 8. The initial potential is NN-only Argonne uig [29j. These results include 
evolved electromagnetic interactions and charge dependences. 
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In previous comparisons between the SRG with flow parameter A and Viowfc 
interactions with momentum cutoff A, the characteristics of the potential and 
its predictions at low energy were found to be strikingly similar in two-body 
calculations when A ~ A. The change in convergence with different A's is 
similar to that shown for different A's in Section [2l Ground-state energies for 
^H, ^H, ^He, and ^Li are plotted against the harmonic-oscillator parameter hfl 
for a range of A^max values in Figs. [9] through [T2l Each flgure shows results for a 
sharp and two smooth Viowfc potentials with A = 2fm~^ and an SRG potential 
with A = 2fm~^. The initial potential is the 500 MeV N^LO potential from 
Ref. [13]. 

We flnd that Viowfc potentials with smoother cutoffs (but the same A = 2 fm^^) 
converge more uniformly and that lower A^max estimates are superior. The lat- 
ter advantage seems to largely disappear as the space gets larger (e.g., see ^Li). 
However, extrapolations may be more robust with the SRG and the smoother 
Viowfc potentials. Note that the SRG is a smooth cutoff low-momentum poten- 
tial. With Td deflned with a single power of T, its characteristics are closest to 
those of a Viowfc potential with exponential cutoff exp[— (A;^/A^)'^] with A ~ A. 

The UCOM framework is an alternative approach to using unitary transfor- 
mations to soften NN potentials and reduce correlations. We can compare the 
observed SRG convergence to that found for the UCOM potential in Ref. [35], 
where results for and ^He starting from AV18 are shown in their Fig. 8. For 
at A^tnax = 6 the UCOM result is about 2 MeV from the converged value 
and at A^max = 12 the discrepancy is about 0.75 MeV. The SRG calculation 
with A = 2fm~^ is 0.5 MeV from converged at A^ax = 6 and only 25keV 
away at A'max = 12. With A = 1.5 fm~^, the discrepancies are 90 and 4keV, 
while with A = 3.0 fm~^, they are 1.9 MeV and 200 keV. If we consider ^He, 
the UCOM result at A'max = 6 is about 2.5 MeV and at A'max = 10 about 
IMeV from converged. For the SRG, the A = 1.5 fm~^ discrepancies are 30 
and IkeV and at A = 3.0 fm~^ they are about 2.6 and 0.4 MeV. Therefore we 
conclude that, in terms of convergence, the UCOM potential from Ref. [32] 
and the SRG A = 3.0 fm~^ potential are roughly equivalent. 



4 Running of Ground-State Energies 

The SRG flow equation guarantees a unitary transformed Hamiltonian and 
exactly the same values for observables only if the entire Hamiltonian is kept. 
However, the commutators in Eq. ([2]) generate many-body interactions as A is 
lowered (e.g., substitute second-quantized operators) and therefore in practice 
there will always be a truncation. But this is also true of the initial EFT 
Hamiltonian, which will already have many-body interactions to all orders 
that evolve with A. However, there is a chiral EFT power-counting hierarchy 
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Fig. 9. Ground-state energy of the deuteron as a function of HQ for a sharp and two 
smooth Viowk potentials with A = 2fm~^ and an SRG potential with A = 2fm^^. 
The initial potential is the 500 MeV N^LO NN-only potential from Ref. [ISj. The 
legend from Fig. [T] applies here. 
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Fig. 10. Ground-state energy of the triton as a function of hQ for a sharp and two 
smooth V^owfc potentials with A = 2fm~^ and an SRG potential with A = 2fm~^. 
The initial potential is the 500 MeV N^LO NN-only potential from Ref. [13]. The 
legend from Fig. [1] applies here. 
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Fig. 11. Ground-state energy of ^He as a function of HQ for a sharp and two smooth 
Mowfe potentials with A = 2 fm^^ and an SRG potential with A = 2 fm^^. The initial 
potential is the 500 MeV N^LO NN-only potential from Ref. [13]. The legend from 
Fig. [1] applies here. 
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Fig. 12. Ground-state energy of Li as a function of h^l for a sharp and two smooth 
Mowfc potentials with A = 2fm~^ and an SRG potential with A = 2fm~^. The initial 
potential is the 500 MeV N^LO NN-only potential from Ref. [13j . The legend from 
Fig. [1] applies here. 



14 



that establishes natural sizes for the contributions from many-body forces. 
This enables truncation at few-body interactions with a controlled error (at 
N^LO this includes a long-range four-body interaction). Therefore the key to 
the usefulness of the SRG for nuclear structure is maintaining this hierarchy 
with a comparable truncation error. In the present calculations, keeping only 
two-body interactions means that the transformations are only approximately 
unitary for 3 and few-body binding energies will change ("run") with A. 
By tracking this running, we can determine the expected net contributions 
from many-body forces in a more complete calculation. The useful range of 
A is where the variations are of natural size (i.e., the same order as the EFT 
truncation errors). 

For the lower A's in the lighter nuclei, our predictions for ground-state energies 
are fully converged. However, in other cases we need to extrapolate the energies 
to A^max = oo. If a labels the hfl values and i the A^max values for each a, a 
possible model for ground-state energies is 

E^, = E^ + A^e-'-''^ , (4) 

where and ba are constants. The goal of a fit to calculations such as those 
in Section [2] is to determine the common parameter Eoo , which is the estimate 
for the ground-state energy extrapolated to iVmax = oo. Presumably this ap- 
proach will work better with larger nuclei for which A^max may become a good 
logarithmic measure of the number of states. 

Given the model in Eq. (jlj), one way to proceed is to minimize the function 
/(Eoo, {6J) = J2iEa^ -E^- A^e-'-''^f/al, , (5) 

where we have allowed the possibility for different weights depending on Ml 
and iVjnax- As posed this is a nonlinear least-squares minimization problem 
with many parameters, for which robust solutions are difficult to find. An 
alternative approach is to recast the problem so that all parameters except 
£"00 are treated as linear parameters. This is possible because our problem is 
strictly variational, so E^i — E^ > for all a and i, which means that the 
logarithm of Eq. (jlj) is well defined. So one considers instead 

log{E^i - E^) = log Aa, - ho,Ni = - h^Ni . (6) 

For each value of E^, this is a linear least-squares problem in the {a^} and 
{ha} parameters. Thus we have a one- dimensional constrained minimization 
problem with the function 

g{E^) = J2{log{Ea, - Eoo) - a„ - 6aA^i)Vf^ai , (7) 



15 



-|- NCSM energ 


ies 


Fit to N = 


2,4,6 


max 


Fit to N - 


2,4,6,8 


max 




max 



Fig. 13. Fits to the NCSM ground-state energies for '^Li at A = 2 fm~^ with HQ 
= 20MeV according to the model in Eq. using the procedure described in the 
text. The dashed hne uses the NCSM energies (pluses) up to A'^^ax = 6, with the 
extrapolated value of E^o shown as the dot-dash line, while the solid line uses up 
to N^s.x = 8, with the dotted line the extrapolated value of E'oo- 



where the {a^} and {b^} are determined directly within the function g by in- 
voking a constrained linear least-squares minimization routine. The constraint 
is the bound E^o ^ min{{Eai}) , where E^o < and "min" means "most neg- 
ative". Again, we allow for weights depending on A^max and/or HQ. 

In the present investigation, we apply Eq. ([7]) with only the HQ value that 
yields the lowest energy in the largest space, weighting different Amax by the 
slope of the energy vs. A^max, and using the spread of results from neighboring 
Hfl values to determine a conservative confidence interval for the extrapolation. 
This procedure is applied for all results in this section. An example of the fit 
to NCSM data is given in Fig. [131 A more complete and systematic study of 
extrapolation in the NCSM will be given in Ref. |33j . 

The extrapolated results for ground-state energies (in MeV) are summarized 
in Table [1] along with the value of hfl (in MeV) for which the energy is minimal 
in the largest spaces. The uncertainties represent conservative confidence in- 
tervals. When not given, the result should be accurate to the digits displayed. 
The and ^He results are fully converged for A ^ 2fm~^. For larger A values 
the energies are evaluated in a large basis (Amax = 48) to ensure con- 
vergence for all A and the "^He results are extrapolated. All of the results for 
A > 4 are extrapolated to some degree, with increasingly small extrapolations 
with decreasing A (and energies for the bare potentials could not be reliably 
extrapolated in the spaces used). 
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Table 1 

Ground-state energies and optimal h^l (in MeV) of light nuclei for SRG-evolved 
potentials. When an uncertainty is given it represents a confidence interval derived 
from an extrapolation to A'^max = oo. Otherwise the energy is converged to the digits 
shown. 
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^He 




°He 




fir ■ 




7t ■ 

'Li 
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hQ 




HQ 


Egs 




Egs 


HQ 


Egs 


HQ 


Egs 


oo 




-7.85 


42 


-26.1(8) 














3.0 


28 


-8.29 


34 


-27.5(3) 


28 


-28(1) 


28 


-31.5(8) 


24 


-38.7(30) 


2.5 


24 


-8.41 


28 


-28.2(2) 


24 


-28.9(3) 


24 


-32.1(3) 


24 


-38.7(20) 


2.25 


22 


-8.47 


24 


-28.6(1) 


22 


-29.4(2) 


22 


-32.5(2) 


22 


-40.3(10) 


2.0 


18 


-8.53 


24 


-28.90 


20 


-30.0(1) 


20 


-33.1(1) 


20 


-41.2(5) 


1.75 


16 


-8.55 


20 


-29.13 


16 


-30.6 


18 


-33.6 


18 


-41.7(4) 


1.5 


12 


-8.48 


18 


-28.86 


14 


-30.7 


16 


-33.7 


16 


-42.0(3) 


1.25 


10 


-8.21 


14 


-27.58 


12 


-29.9 


12 


-32.9 


12 


-41.1(2) 


1.0 


8 


-7.63 


14 


-24.80 


10 


-27.4 


10 


-30.4 


12 


-37.8(2) 



If we calculate the energy of a nucleus using the SRG-evolved NN interaction 
only, the running of the energy is the running of the net three-body (and 
higher-body) contribution to the energy. Over a wide range of cutoffs, many- 
body interactions generated in the RG will be the same scale as the starting 
3N interactions in the chiral EFT. Figures [T^ to [TBI show the NN-only running 
for ^H, ''He, ^He, and ^Li. The energy plotted is either the converged value or 
extrapolated from the best HQ with an error bar (i.e., a confidence interval) as 
described above. Figures [T71 to [20] are binding-energy correlation plots between 
various nuclei. 

The pattern of SRG running for ^He, ^He, ^Li, and ''Li is qualitatively similar 
to the pattern for shown in Fig. [HI That is, a slow increase in the binding 
energy as A decreases to A = 1.5-1.8 fm~^ and then a steep decrease. The two 
SRG curves in Fig. [T^ show the difference between evolving the neutron-proton 
interaction only and evolving the full NN interaction. The SRG pattern is also 
qualitatively similar to the running for Viowfc potentials [8|[5] . As observed for 
^owfc [m], the size of three-body contributions is natural according to chiral 
EFT power counting. This is true for the net energy within the entire range of A 
shown in the figures, but future calculations are needed to investigate whether 
the individual contributions from long-range and short-range 3N forces remain 
natural. The smoothness of the running adds credibility to the extrapolations 
and the associated estimates of confidence intervals; when the estimate has 
been improved (e.g., from a larger A^max result), the running and correlation 
plots have become more regular. 
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Fig. 14. Plot of the ground-state energy of the triton vs. A (if SRG) or A (if Viowfc) for 
potentials evolved from the 500 MeV N^LO NN-only potential from Ref. [13]. The 
arrow marks the experimental binding. The two SRG curves show the difference be- 
tween evolving the full NN interaction and evolving the neutron-proton interaction 
only. 
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Fig. 15. Plot of the ground-state energy of ^He and ^He vs. A for potentials evolved 
by the SRG from the 500 MeV N^LO NN-only potential from Ref. Conservative 
error bars have been included with the larger A's, for which an extrapolation is 
needed. The arrow marks the experimental binding. 
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Fig. 16. Plot of the ground-state energy of ^Li and ^Li vs. A for SRG potentials 
evolved from the 500 MeV N^LO NN-only potential from Ref. [13]. Conservative 
error bars have been included with the larger A's, for which an extrapolation is 
needed. The arrow marks the experimental binding. 
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Fig. 17. Tjon line for SRG potentials evolved from the 500 MeV N^LO NN-only 
potential from Ref. [13j (with A in fm~^). Conservative error bars have been included 
with the larger A's, for which an extrapolation is needed. The result for ^He with 
the bare potential is from Ref. [36] . The line is a fit to NN potentials from Ref. [37] . 
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Fig. 18. Binding energy correlation plot between ^He and '^H for SRG potentials 
evolved from the 500 MeV N^LO NN-only potential from Ref. [131 (with A in fm"^). 
Conservative error bars have been included with the larger A's, for which an extrap- 
olation is needed. 




E^C^H) [MeV] 

Fig. 19. Binding energy correlation plot between ^Li and for SRG potentials 
evolved from the 500 MeV N^LO NN-only potential from Ref. p3] (with A in fm"^). 
Conservative error bars have been included with the larger A's, for which an extrap- 
olation is needed. 



20 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



34 



A,=1.5 




33 



A,=1.25 




A,=2.5 



i 



h=2.25 



32 



I — n — [ 



Expt. 



A,=3.0 



31 



X X=1.0 




Ej^(^He) [MeV] 



Fig. 20. Binding energy correlation plot between Li and He for SRG potentials 
evolved from the 500 MeV N^LO NN-only potential from Ref. [IS^ (with A in fm~^). 
Conservative error bars have been included with the larger A's, for which an extrap- 
olation is needed. 

The correlation plot for ^He vs. shows the expected movement along the 
Tjon line, with a slight loop closing to the left as A decreases. The latter be- 
havior is amplified in the ^He vs. and ^Li vs. plots. The plot of ^Li vs. 
^He is very close to linear. Tjon lines are also found in coupled cluster calcu- 
lations of ^^O, ^^O, ^''O, ^^N, and ^''F [38]. From these plots one sees some A's 
for which predicted ground-state energies are quite close to experiment (e.g., 
A = 2.25 fm"^). These may be good choices for calculations for which small 
net three-body contributions are desirable. However, we emphasize that this is 
not a global trend. Indeed, nuclear matter does not saturate with only NN in- 
teractions at such A's (except perhaps at very large density). Thus calculations 
in heavier nuclei (e.g., ^^O, ^"Ca, . . . ) will be increasingly overbound. 



5 Other Observables 

In this section, we show some limited results of excitation energies and radii 
for ^Li. Because we expect sensitivity to three-body contributions, we include 
these results primarily to assess convergence trends as a baseline for future 3N 
calculations and to compare to previous results from Ref. [TS]. 



21 




[MeV] 

Fig. 21. Excitation energies of the lowest natural-partiy states of '''Li as a function of 
hn for A = 2.0 fm"^ The initial potential is the 500 MeV N^LO NN-only potential 
from Ref. [13j. The horizontal dotted lines are the experimental values while the 
vertical dotted lines mark the optimal hi} value for the ground-state energy (middle) 
and the range for which estimates are close to this. 




[MeV] 

Fig. 22. Excitation energies of the lowest natural-parity states of ''Li as a function of 
hn for A = 1.5fm"\ The initial potential is the 500 MeV N^LO NN-only potential 
from Ref. [l3j. The horizontal dotted lines are the experimental values while the 
vertical dotted lines mark the optimal hil value for the ground-state energy (middle) 
and the range for which estimates are close to this. 
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The spectrum of low-lying natural-parity excited states in Li is shown in 
Figs. EI] and [22] for A = 2.0 fm"^ and 1.5 fni'^ The excitation energy is plotted 
against hQ for different A^max values to facilitate comparison to similar plots 
from NCSM Lee-Suzuki calculations based on the same potential [TS]. The 
vertical dotted lines mark the range of HQ that has the lowest ground-state 
energy estimates, with the middle line marking the best energy estimate for 
the ground state. The higher A/'max predictions are flatter as a function of HQ, 
as expected, and A = 1.5 fm~^ converges faster than A = 2.0 fm~^. However, we 
observe that the curves are not as flat as in the corresponding plot in Ref. [18j. 
The spread from A^max = 2 to A'max = 6 is not large, but also appears larger 
than in Ref. [H]. 

In Figs. [2S] and [M] we plot the point proton radii of ^He and ^Li as a function 
of HQ, with the different A^max curves having the same legends as in Figs. [T] 
and O These plots are made using the bare operator for at all A, rather 
than an evolved operator. (Note that for the proton is not an observable.) 
Since samples primarily long distances, its matrix elements depend weakly 
on the cutoff [5]. For the deuteron, the rms radius changes from the X = oo 
value by less than 1% by A = 1.5 fm~^ and by about 4% by A = 1.0 fm~^ 
Our observations are consistent with the investigations reported recently on 
the range- dependence of another renormalization scheme [39] . 

For ^He, we observe that with decreasing A the curves for a given A^max get 
flatter; for A = 1.5 fm~^ our A^max = 12 results are essentially independent of 
hfl. The different A^max curves cross each other in a decreasing range of HQ as 
A decreases. For A = 1.5 fm~ , all of the curves for both nuclei intersect at the 
same point. If we consider the intersection of the three highest A^max results, 
it is close to (but generally slightly below) the HQ for the lowest ground- 
state energy estimate. A similar observation has been reported for the NCSM 
Lee-Suzuki calculations in Ref. [TS]. The trend is consistent with the idea 
that if you have the optimal HQ for a given nucleus (which scales with A~^^^ 
and inversely with the average radius squared), then the radius should be 
minimally sensitive to changes in the basis size. For ^Li, the intersection point 
is well below the optimal HQ for the energy except at very low cutoffs; the 
radii for larger A are not converged in the largest spaces considered here. 

In Fig. [211 the dotted lines indicate the radii from the NCSM Lee-Suzuki cal- 
culations of the bare N^LO potential without 3N interactions (middle) and 
with two different 3N fits (top and bottom) [IB]. The spread of these predic- 
tions can be taken as an estimate of the expected 3N contribution for various 
A's. The hQ dependence here is steep, but if we take the intersection point to 
determine the NN radius, then the results for A = 3,2, 1.5 fm~^ are all com- 
patible with the radii from Ref. [I8]. Moreover, the trend in the radius with 
decreasing A parallels the running of the energy with A (see Fig. [TB|1 . 



23 



B 1.7 



B 1.7 




o 

Oh 



15 20 25 
bQ. [MeV] 



25 30 
iiQ [MeV] 

Fig. 23. Point proton radii of ^He as a function of HQ at four different values of A 
(oo, 3, 2, 1.5 fm"^). The initial potential is the 500 MeV N^LO NN-only potential 
from Ref. [iSj. The legend from Fig. [T] applies here. 
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Fig. 24. Point proton radii of ^Li as a function of HQ at four different values of A (3, 
2, 1.5, Ifm"^). The initial potential is the 500 MeV N^LO NN-only potential from 



Ref. [13]. The three dotted lines are radii from NCSM Lee-Suzuki calculations using 
the bare N^LO potential in larger spaces pS]. The middle (black) line is without 
3N interactions and the others are with two different 3N force fits [H]. The legend 
from Fig. [5] applies here. 
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6 Summary 

We have performed NCSM calculations that retain the variational principle 
using RG-evolved potentials in momentum representation, focusing on the 
SRG approach that has recently been applied to nuclear structure physics. 
While we have not performed an exhaustive test of alternative potentials, we 
have found no sign of any special dependence on the initial potential. The 
convergence of ground-state energies improves rapidly with decreasing A and 
with no unnatural contributions from omitted three-body forces for these light 
nuclei. Furthermore, all methods of evolution to low momentum give similar 
convergence improvements, with smoother behavior associated with smoother 
cutoffs. 

With these calculations we have set the stage for including 3N interactions |23j . 
Just as varying A in Figs. to [T6] shows the scale of omitted many-body con- 
tributions, similar figures when 3N is included will quantify the net impact 
of higher-order interactions. In the short term, three-body forces will be in- 
cluded in the form of N^LO chiral 3N interactions, with the parameters fitted 
(see Ref. [8] for motivation). The development of evolved 3N potentials is pro- 
ceeding in parallel and should be available in the near future. Once we have 
consistent 3N interactions, we will be able to establish lower bounds to the- 
oretical error bars for light nuclei by calculating with different N^LO cutoffs 
and assess the impact of uncertainties in the input chiral 3N interactions. 
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